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We present two methods to include the asymptotic domain of a background spacetime in null di- 
rections for numerical solutions of evolution equations so that both the radiation extraction problem 
and the outer boundary problem are solved. The first method is based on the geometric conformal 
approach, the second is a coordinate based approach. We apply these methods to the case of a 
massless scalar wave equation on a Kerr spacetime. Our methods are designed to allow existing 
\ codes to reach the radiative zone by including future null infinity in the computational domain with 

. relatively minor modifications. We demonstrate the flexibility of the methods by considering both 

' Boyer-Lindquist and ingoing Kerr coordinates near the black hole. We also confirm numerically 

fSJ I predictions concerning tail decay rates for scalar fields at null infinity in Kerr spacetime due to Hod 

for the first time. 
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I. INTRODUCTION 



In numerical calculations of radiative fields, it is common to artificially truncate the computational domain by 
I , introducing an outer boundary into the spacetime. It is well known that this practice introduces conceptual and 
^jPj' operational problems, specifically the outer boundary and the radiation extraction problems. [U, S S S S S B S 
[Tol . [ill . [T^ . ITsl . [3, [l^, 13 [l3l ■ Each incremental development has taken considerably more effort, while a clean solution 
to the problems is available on a geometrical level, namely including the asymptotic domain on the computational 

grid di, m, M M m, m, m . 



• First attempts on including the asymptotic domain on the numerical level have been based on using a compactifying 
' coordinate along outgoing null surfaces [25l . [26j . This allows one to include null infinity in the computational domain 
. where no boundary conditions are required and radiation is naturally and unambiguously measured by idealized 
' observers. In addition, the numerical evolution is very efficient as the grid follows outgoing radiation. Setting the 
. numerical outer boundary to null infinity solves the aforementioned problems in a numerically efficient manner. 
' Unfortunately null foliations can be inconvenient, especially because they tend to develop caustics in dynamical 
^\ [ spacetimes [27| . For solving the outer boundary and radiation extraction problems, however, we only need access 

■ to null infinity, whereas the interior foliation can be chosen arbitrarily. This observation suggests that one could 

■ match an arbitrary spacelike slice in the interior to an asyrnptotic domain that extends to null infinity. This has 
first been tried using Cauchy-characteristic matching [H, [H, [2^, [sOl ■ There, one uses an evolution based on Cauchy- 
type foliations in the interior and communicates the in- and outgoing characteristic variables along a timelike surface 
with a characteristic evolution in the asymptotic domain. This method introduces numerical difficulties due to the 
presence of a nonsmooth matching along which the causal nature of the slicing changes from spacclike to null. While 
Cauchy-charactcristic matching is a promising approach that requires further study, in this paper we will explore an 
alternative approach to the problem. 

Instead of matching spacelike slices to null slices along a timelike surface, here we use everywhere spacelike slices 
that approach null infinity in the asymptotic domain. The latter are called hyperboloidal surfaces as their behavior 
is similar to the standard hyperboloids in Minkowski spacetime [l9| . They are more flexible than null surfaces 
because the only local condition restricting their choice is that they are spacelike. The asymptotic condition that 
they approach null infinity does not restrict the type of spacelike surfaces in the interior. It is favorable for numerical 
applications within the hyperboloidal approach that a g au g e is chosen in which the location of future null infinity 
(scri) is independent of time, i.e. scri fixing [3l|, [H, [U, [SJ, [s^, [36l . Numerical experiments with scri fixing have 
already been made successfully in spherically symmetric spacetimes |37l . Issl . Isol . liol . l4ll . 

In this paper we discuss the case of a fixed Kerr black hole background from the point of view of the hyperboloidal 
approach with scri fixing. One might think that the characteristic approach should be sufficient for calculating radiative 
fields at null infinity when a background has been given. Indeed, for given Minkowski and Schwarzschild backgrounds, 
null coordinates are very useful and have been regularly applied in numerical calculations [26. . ,43, , ,441 , |45| | . In Kerr 
spacetime, however, they are difficult to deal with. There has been an ongoing effort to find a metric representing 
the Kerr geometry in Bondi-Sachs form that can be used in numerical computations [1^, H^, 113, E^, , but to our 
knowledge no numerical calculations including null infinity could be presented so far. 

The difficulties with the characteristic approach in Kerr spacetime arise, in part, from the fact that null coordinates 
are very rigid, i.e. the coordinate freedom for choosing a null surface is more restricted than for a spacelike surface. 
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The flexibihty of hyperboloidal surfaces, on the other hand, allows us to evolve radiative fields on a Kerr spacetime up 
to and including null infinity while keeping, for example, standard Boyer-Lindquist or ingoing Kerr coordinates near 
the black hole ^. We have chosen these foliations as they are the most common ones and have a differing qualitative 
behavior near the horizon. In principle, any foliation near the black hole can be chosen. The construction of such 
hyperboloidal surfaces in Kerr spacetime is relatively new (36j . In this paper we show that it is also amenable to 
numerical calculations. 

We discuss two methods for the numerical implementation of the hyperboloidal matching idea on the example of 
a scalar wave equation. In both methods matching includes certain choices of rescaling, time transformation and 
coordinate compactification. The first method is the well-known conformal method as introduced by Penrose [12, [2^ . 
It requires knowledge of the conformal transformation behavior of the equation one is trying to solve. The second 
method is an essentially equivalent method that may be more appealing to researchers unfamiliar with conformal 
techniques. It considers the equation to be given in a certain coordinate system and applies the above-mentioned 
transformations (4^ . These methods allow us to solve the equations of interest including null infinity and can be 
applied within the hyperboloidal approach independent of matching. 

We test our imp lementations by reproducing previously obtained quasinormal mode frequencies [5ll . [5^ and tail 
decay rates [13, Ej, HI] • Our numerical tests focus on scalar fields but our method of including null infinity can be used 
with a wide variety of systems of evolution equations on a Kerr, Schwarzschild or Minkowski background that admit 
asymptotically flat solutions. In particular, the application of our techniques to calculate other types of radiative 
fields such as gravitational perturbations should be straightforward. 

The paper is organized as follows. In Sec. II we present the main idea of hyperboloidal matching that allows us 
to include null infinity in the computational domain while using standard coordinates near the black hole. Sec. Ill 
presents the conformal method and applies it using ingoing Kerr coordinates. Sec. IV presents the coordinate based 
method applied to Boyer-Lindquist coordinates. The resulting equations are then solved in Sec. V using standard 
numerical techniques. Beyond confirming well-known predictions with our code, we also study decay rates for scalar 
perturbations of a Kerr black hole near and at null infinity confirming predictions due to Hod [sH] for the first time. 
We conclude with a short discussion of our results. 



II. THE MAIN IDEA 




As mentioned in the introduction, the construction of null coordinates in a Kerr background is relatively complicated 
J4^ whereas hyperboloidal foliations that are useful for numerical purposes can be given explicitly [36j . 
36[ , we construct hyperboloidal foliations that coincide in an interior domain exactly with ingoing Kerr or 
Boyer-Lindquist coordinates. 

The ingoing Kerr and Boyer-Lindquist representations of the Kerr metric will be given in terms of coordinates that 
we denote by (t, r, y, (p). Here, t is the time coordinate of the corresponding metric, r is a radial coordinate along level 
sets of t and y is an angular coordinate defined by cos'd where and ip are the angular coordinates on a sphere. In 
order to match a spacelikc surface to null infinity in a stationary spacetime, we introduce a new time coordinate r of 
the form 

T = t-h{r). (1) 

The function h{r) is called the height function. The advantage of the above transformation is that it keeps the time 
direction invariant. The stationary Killing field that is timelike outside the event horizon has the same coordinate 
representation in both t and t for any choice of h(r). It is given by dt or, equivalently, by dr- This implies that 
the Kerr metric in the new coordinates is manifestly stationary and natural observers that follow integral curves of 
the stationary Killing field are again given by surfaces of constant spatial coordinates. A natural consequence of the 
choice ([1]) is that null infinity is at a fixed spatial coordinate location with respect to a time- independent compactifying 
coordinate a long level sets of r. This considerably simplifies the numerical implementation as well as the analysis of 
the solution [saEI- 

A compactifying coordinate can be introduced by 

r=^, with n = n{p). (2) 



^ We use the naming convention for the coordinatization of Kerr spacetime as in [soll . 
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This choice imphes that the zero set of corresponds to infinity in terms of the physical coordinate r. The choice 
of the conformal factor can be made in such a way as to ensure that future null infinity (denoted by is at the 
coordinate location p = S, for example, by setting f2oxt(p) = 1 — p/S. In an interior domain near the black hole we 
want to use the standard physical coordinates, which implies flintip) ~ 1- There is a transition between the interior 
and the exterior domain where we match the two functions. Many choices are possible on the transition domain. Here 
we set 

{^^int(/3) for p<Pint, 

nintip) e-^P-P"^^^"/^P-P^-^^" + (l - e-(''-'''"«)'/(P-PexO=) for p,^^, < p < p^^^, (3) 

^^oxt(p) for P>Pcxt- 

Hyperboloidal foliations that coincide with standard Boyer-Lindquist or ingoing Kerr foliations in the interior can 
be constructed by suitably choosing the derivative of the height function dh/dr =: H. Here we choose 

{0 for p < Pint, 

-ffcxt(p) (l - e-(''-'''-)'/(''-^-')') for Pint <P< Pext, (4) 

Hcxt{p) for P > Pcxt- 

The choice H = along with = 1 near the black hole ensures that the foliation and the coordinates are not changed. 
The function Hext{p) is different for Boyer-Lindquist and ingoing Kerr foliations [36j . Its choice will be discussed in 
later sections. 



III. THE CONFORMAL METHOD 

The conformal method to study the asymptotic behavior of radiative fields was introduced by Penrose [13, [H, [s^ . 
The example of the scalar wave equation was the first system where the hyperboloidal initial value problem was studied 
(s^. |58|. In this section we will present the conformal method for the scalar wave equation on a Kerr spacetime. 



A. Conformal transformation of the scalar wave equation 

Our notation follows [4l|, where the conformal method was applied to calculate tail decay rates for scalar and 
Yang-Mills fields at null infinity in a Schwarzschild spacetime. 

Here we are interested in solutions to the scalar wave equation Dcp — 0, for a scalar field in a Kerr spacetime 
with metric g. The coordinate transformations ([T]) and ^ result in a representation of the Kerr metric g that is 
singular at fl — 0. To include null infinity in a regular way we need to rescale the Kerr metric g with Q,^. The Kerr 
spacetime is weakly asymptotically simple, implying that with a suitable choice of the transformations ([1]) and ([2|) the 
rescaled metric g = fl^g is smoothly extendable through null infinity where we have {fl = 0, d^l ^ 0}. The following 
conformal transformation rule holds 

n - = (^fi - ii?^ ^, with 0=^, (5) 

where R and i? are the Ricci scalars of the rescaled and the physical metrics g and g respectively. The scalar wave 
equation with respect to the rescaled metric then becomes 

□0 - = 0. (6) 
Note that the above equation is invariant under both coordinate and conformal transformations. 

B. First order reduction 

The wave equation ([6]) can be written in any coordinate system {a;''} as 

= Y^d^<\> + where g^^Vl^. 



4 



To bring it to first order form we introduce the auxiliary variables 4>^ := d^cj). Assuming that the coordinate 
denotes the time direction we get the following system of evolution equations: 

do(f>a = 9a 00, (7) 
9o0 = 00- 

The above system is symmetric hyperbolic [s^. Its coefficients are calculated with respect to the conformal metric 
g written in coordinates {t, p,y,ip}. In the case of a Kerr spacctimc we can apply the separation of variables 
(j){T, p,y,ip) = (f){T^ p^y)e~'^^^ and study solutions for each k mode independently, because azimuthal modes decouple 
in the presence of axisymmetry. Here, we will choose initial data corresponding to a single azimuthal mode with 
/c = 0. In this case, the indices a and b span over p and y. 



C. Hyperboloidal compactification in ingoing Kerr coordinates 



In the following, we present a hyperboloidal compactification of the ingoing Kerr representation of the Kerr metric. 
The ingoing Kerr slicing is one of the most commonly used slicing conditions for numerical calculations on Kerr 
spacetime. Its level sets allow us to apply an excision technique inside the black hole horizon. Asymptotically they 
approach spatial infinity. For their relation to the Boyer-Lindquist coordinates see (50l |. Following usual notation we 
introduce 

S:=r2+a2y^ A := r'^ + - 2mr. (8) 

where m and a are the ADM mass and the angular momentum of the Kerr spacetime, respectively. After transforming 
to compactifying coordinates we will use S := fi^E and A :— fi^A. 
The Kerr metric in ingoing Kerr coordinates reads 

.giK = - ( 1 ^\dr + ^-dtdr ^{1 - y^) dtd^p - 2a{l - y^) \A + \ drdip + 



j„2 , ^ J 2 , f^2 . 2 , 2ma'^r{l-y^) 



1 + ^] dr^ + dy' +[r' + a' + ^ (1 - y') dip' . (9) 



This metric gives the following simple evolution equation for dtcfit- 

(S + 2mr) dt<j)t = 4mrdr(f)t + A c^^^r + (1 - y'^)dy<f)y + 2m(f)t + 2{r ~ m) (f)r - 2y4>y. (10) 
Applying the coordinate transformations H]) and ^ to the above metric and rescaling with il^ gives 

5iK = -n^(l-'=^)dr^-2L (H-'=^il + H)) drdp-'-^il-y^)drd^ + g,,dp^~ 



2La{l ~y')(l + ^(1 + H)] dpd<p + -^dy' + ( p' + a'n' + ^"^"'^^'(l ' ) (i _ y2)^^2 _ ^^^^ 



Here, L := — pdp^l, and 



L f „2 2mpQ. 



The radial component of the spacetime metric appears to be singular, but with a suitable choice for the derivative of 
the height function it acquires an explicitly regular form. Following [36| we make that choice in the exterior domain 
to be 

g„..i + l;!!g + («'"'-fg<)"' , ,12) 

p p^ 

where CiK is a free parameter. Only the first two terms are required for the regularity of gpp. The third term (involving 
the free parameter Qk) is taken for numerical purposes as will be discussed later. The appearance of such a free 
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parameter for hypcrboloidal foliations should be contrasted with asymptotic solutions to eikonal equations used for 
asymptotic constructions of null hypersurfaces, where no such parameters are allowed [46| . 

With the above choice for the derivative of the height function the metric is manifestly regular at J^+. By setting 
f^cxt = 1-/5/5" we get 

giK\j'+ ^ -2dTdp + 2-^ dp^ - 2a(l - 2/2) dp dip + dy^ + (1 - y^) dip^. 

The coordinate speed of outgoing characteristics at null infinity reads c+ = —2grp/gpp', for the above metric this is 
c+ = S'^/C'^^. We see that the coordinate speed of characteristics depends on the coordinate location of null infinity 
and an additional free parameter CjK- The latter is related to the asymptotic value of the mean extrinsic curvature of 
our surfaces. This is the reason why we included the third term in the expansion (|12p : the parameter Cik is important 
if we want the coordinate location of to be large while avoiding strong restrictions on the allowed time step due 
to a large value of the coordinate speed of outgoing characteristics. 

IV. THE COORDINATE BASED APPROACH 

We can recognize essentially three steps in the hypcrboloidal compactification of Kerr spacetime: a time trans- 
formation ID), a spatial compactification and a rescaling In the previous section we calculated a conformal 
metric that has a smooth extension through null infinity and studied the conformal transformation behavior of the 
underlying covariant differential equation. In this section we will instead go through the above three steps at the level 
of the partial differential equations written in some coordinate system. We will follow this approach in an explicitly 
noncovariant notation. The rationale behind this is to have an alternative to the geometric approach, which might 
be more straightforward to apply for those not familiar with conformal techniques. This approach may also be useful 
in cases where the conformal transformation behavior of the equation of interest is difficult to calculate [4^ . 

A. The transformations 

Suppose that we are given a linear, second order partial differential equation for a function in two spatial 
dimensions 



dt(t>t = A'^drCjjt + A"dr(t)r + dyCjiy + SVt + B^(j}r + By(l)y + C0, (13) 

where (j)p_ := dp(j). Assuming that the coefficients are independent of <, we write the above equation under the following 
transformations: 

(i) Time transformation: We introduce a new time coordinate as in ([T])- Equation then becomes 

(1 - A"H^ + A*''H) dr4>r = (i*'' - 2HA'')dr4>r + A''''dr4>r + dyij)y + 

+ {B' -B''H-A''''drH)4>r + B''4>r+By4>y+C4>, (14) 

(ii) Spatial compactification: We introduce a compactifying radial coordinate as in ([2]). We set 

J := g = ^, J' := ^ = ^i2idpn)L + pno^n), 

with L = il — p dpfl as before. Equation ([T^ becomes 

dt^t = A'^-JOp^t + A^'-J^ dp4>p + Ayydy4>y + B'~4>t + (BV + i'-V) 4>p + By4>y + C4>. (15) 

(iii) Rescaling: We define the rescalcd evolution variable ((> := f{r)(f>. Equation (fT3l) becomes 



dtcj)t = A'^-drC^t + A'-'dAr + Ayydy(jyy + \B' - A'" ^ j 0t + \B'' - 2A''^ ^ j 0, + + 

C-ff^+A-(2^-^-^))^. (16) 
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Combining the three steps above gives us finaUy an evolution equation of the form 

= A'Pdp(l)r + APPdp(j)p + Ayydy(l)y + B^d^r + B^^p + + C(^. (17) 

This method is essentially equivalent to the conformal method. Equation ([7]), which we obtained with the conformal 
method, can also be derived by using the coordinate based approach just presented and making consistent choices 
of rescaling and coordinate transformations. This has been used as a check for the final equations that the have 
numerically implemented. 

Note that in the coordinate based method we do not need to rescale the function (j) with Q.. We can choose to 
rescale it with any function that behaves asymptotically like r, and read off the radiation field directly at null infinity. 
In fact, a common approach in numerical calculations is to evolve the variable (p := rcf). To see how this relates to the 
conformal rescaling of the previous section, make a simple choice for the conformal factor f2 = 1 — p. In terms of the 
physical coordinate r we have have 

P ^ 1 

1 — p 1 + r l + r 

Therefore and due to the fall-off behavior of (/), the rescaling = (1 + r)(j) is asymptotically equivalent to the 
rescaling r0 in first order. 



B. Hyperboloidal compactification of the scalar wave equation in Boyer-Lindquist coordinates 

We have used the above method to transform the scalar wave equation on a Kerr background written in Boyer- 
Lindquist coordinates. The Kerr metric in those coordinates reads 

ffBL = - 1 - ^ di' ^(1 - y') dtdip+^ dr' + 2 dy^ + (1 - y^) d<p2_ 

Instead of writing the conformal Kerr metric we derive the scalar wave equation on this background and work directly 
with that equation. It is common practice to introduce the tortoise coordinate for the numerical implementation of 
the wave equation because the Boyer-Lindquist coordinates do not give a foliation of the event horizon. The tortoise 
coordinate is defined by 

dr^ + o? 

dr A 

The scalar wave equation in Boyer-Lindquist coordinates can be written as 



dt(j)r, = dr,4>t, dt<j)y = dy<j>t, df4> = 4>t, 



where we have defined := (r^ +a^)^ — a^A(l — y^). From the above equation we calculate the equation for a rescaled 
variable (p -.^ r (p following the prescription given by Eq. ([T6|). Then we introduce a time function r as in Eq. (|14p. 
such that the level surfaces of t approach null infinity. The Kerr metric in Boyer-Lindquist coordinates approaches a 
Schwarzschild metric in standard Schwarzschild coordinates; therefore, we may assume that the asymptotic form of 
H needs to be chosen as for the Schwarzschild spacetime in standard Schwarzschild coordinates. Remember, however, 
that the asymptotic Schwarzschild metric in null directions in the tortoise coordinate is of standard Minkowski form. 
Therefore, we can choose for H the expression for the standard hyperboloids in Minkowski spacetime [s^. We set 

^ Vrilci^ ^ T^A^- ^^^^ 

For the second part of the above equation we used r» = p/fl. Here Cbl is a free parameter that determines the 
asymptotic extrinsic curvature of the hyperboloidal surfaces as in the case of Minkowski spacetime [36l . l60j . As in the 
ingoing Kerr case, the value of Cbl needs to be chosen close to the value of the coordinate location of null infinity, S, 
in order to prevent the coordinate speeds of outgoing characteristics from becoming too large. 

For the spatial compactification ^ we choose the function as in the ingoing Kerr case. Here, fl determines the 
spatial coordinate compactification and does not acquire the meaning of a conformal factor for the spacetime metric. 

The above choices fix the transformations and result in an explicitly regular set of equations each of which has the 
form (|17p . Their final form is rather lengthy but straightforward to calculate and will therefore not be presented. 
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V. NUMERICS 

As a test of the two approaches to include nuU infinity in numerical simulations by spacelike matching " standard" 
Cauchy hypersurfaces to hyperboloidal slicings here we implement them on a 2 + f code solving a scalar wave equation 
on a Kerr spacetime. Near the black hole we use ingoing Kerr and Boyer-Lindquist coordinates. 

The conformal method is applied for the hyperboloidal compactification in the ingoing Kerr case, and the coordinate 
based approach in the Boyer-Lindquist case. We note that there is no specific reason for this choice. Other coordinate 
choices are possible as well. The equations we solve numerically are given in Eq. ([7]) with respect to the metric (jlip 
and in Eq. p8|) with the transformations (fT4|) . (fT5| . and (fT6)) . The asymptotic form of the height functions for ingoing 
Kerr and Boyer-Lindquist foliations are and p^ . respectively. 

In both cases we prescribe initial data of the form 

0(O,p,y) = O, 0,(O,p,2/) = e-('^-'-°)'/'^'l^o (20) 

where Yiq denotes the familiar spherical harmonics of degree I and order fc = 0. The data are of compact support. 
Wc set (j)r{0,p,y) = for p > pint- Regarding the parameters in the initial data we choose rg = 5 and a = 2 with 
mass m — 1. 



A. The numerical method 

We use method of lines with fourth order Rungc-Kutta time inte grat ion. For the angular derivatives we use a 
pseudospectral collocation method with Legendre polynomials, as in [5l], which is both efficient and automatically 
takes care of the coordinate singularities at the poles. We use finite differences in the radial direction. The reason 
why we do not use spectral differencing in the radial direction is the nonsmooth matching ^ and We wish to 
remark, however, that we could, in principle, apply spectral differencing by using a multipatch approach in which the 
matching boundaries coincide with some of the patch boundaries. We did not follow this approach only due to the 
simplicity of the finite difference method. 

In the ingoing Kerr case the inner boundary is a spacclikc surface inside the event horizon so that one can use the 
excision technique. The outer boundary is at null infinity. We apply one-sided finite differencing at both boundaries 
using the same order operators at the boundaries as in the interior [6l|. In the Boyer-Lindquist case we set the inner 
boundary very close to the black hole. We apply the outgoing boundary condition (ingoing with respect to the black 
hole) at the inner boundary, which we set at = — 40to. We use fourth order accurate finite differencing operators in 
the Boyer-Lindquist case instead of the eight order operators used in the ingoing Kerr case, as the outgoing boundary 
condition is simpler to apply with respect to fourth order operators. 

We add Kreiss-Oliger type artificial dissipation to the evolution equation for the time derivative of to suppress 
numerical high-frequency waves (6^ . For a 2p — 2 accurate scheme we choose an operator (Q) of order 2p as 

where h is the grid size, D± arc the forward and backward finite differencing operators and e is the dissipation 
parameter which we set as e = 0.03. 

One major source of error is the matching between the inner coordinates and the coordinates of the hyperboloidal 
compactification. The matching is performed using Eqs. ([3]) and ([?]). We found that the matching along pi^t introduces 
large numerical errors in the solution that converge with the order of the finite differencing method. The fact that 
the numerical errors are large near pint is to be expected because the matching is not analytic at pint whereas it is 
analytic at pext. In our simulations we set pint to be at around 20m in both cases. This is a point where improvement 
may lead to a more efficient code. The outer matching introduces no problems. 

We put the numerical outer boundary at p = S" = 100 in both cases. The numerical outer boundary corresponds 
to future null infinity. Its coordinate location effects the coordinate speed of characteristics as discussed before. 
Experiments with different coordinate locations of future null infinity delivered qualitatively similar results. An 
important point with respect to the coordinate location of the outer boundary is related to the choice of parameters 
that determine the hyperboloidal foliation. The free parameters Qk and Cbl need to be both of the order of S as 
discussed in Sees. IIII Cl and HV Bl 

The convergence properties of our codes can be seen in Fig. [TJ For these plots, a three level convergence analysis 
has been performed for an I = initial data using 1600, 3200, and 6400 grid cells in the radial direction and 4 points 
in the angular direction where spectral differencing was applied. We used quadruple precision for the runs because 
the numerical error was below machine accuracy at late times when double precision was used. The time step for the 
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Figure 1: Convergence in the L2 norm in the radial direction for ingoing Kerr (left plot) and Boyer-Lindquist (right plot) 
coordinates near the black hole. The convergence factor Q is calculated by Q = logj !!tmad_t,high\\'^ where the norm is taken 
over the whole grid at a given time. We see that in the ingoing Kerr case the convergence factor attains its value 8 after a 
transient phase because numerical errors in the time integration dominate as observed in [43 |. In the Boyer-Lindquist case the 
convergence factor is 4. The irregularities in the convergence factors are due to matching. 



ingoing Kerr case was dt — 0.05 while for the Boyer-Lindquist case it was dt = 0.1. The reason why we can choose a 
larger time step in the Boyer-Lindquist case is because the coordinate speeds of characteristics are lower. 

We see that, in the ingoing Kerr case, the convergence factor settles down at around 8 after an initial transient 
phase in which it starts at 4, the order of the time integration. The reason is that in this initial phase the error is 
dominated by the time integration. This is the same behavior as in To get a convergence factor of 8 from the 
beginning a much lower time step needs to be chosen which is not practical for late-time tail calculations [i^l • In the 
Boyer-Lindquist case where both the time integration and the spatial discretization is of fourth order the convergence 
factor is 4. The irregularities in the plots are due to matching. When a wave package passes through the matching 
boundary pi^t, the error increases significantly but in a convergent fashion. 
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Figure 2; Quasinormal mode ringing of Kerr spacetime with m = 1 and a — 0.5 excited by a Gaussian initial data with I = 2 
mode as measured by an observer located along r = 3m (the solid red curve) and J^"*" (the dashed blue curve) in ingoing Kerr 
coordinates. 

As a next test we calculate a fundamental quasinormal mode frequency. The quasinormal mode ringing for initial 
data with I — 2 can be seen in Fig. [2l We plot here the ringing as measured by an observer at r = 3m and at J^"*" in 
ingoing Kerr coordinates. We observe that the duration of ringing along is shorter. Therefore we use the signal 
measured closer to the black hole for the fitting. The frequencies are obtained by fitting the signal to the formula 

0(r) = Ae-'^^^sin(LJiT-hv3). (21) 

Here, ui and W2 ai'c the mode frequencies, A is the amplitude and (p is the phase of the wave signal. The fitting 
is performed using a simple least squares method on the interval r € [60m, 150m]. The error in frequency is rather 
dominated by the fit than by numerics. We find uJi = 0.491971 and uj2 = 0.094653. These numerical values are 
very close to those obtained by using Leaver's continued fraction method [H, [63, which read loi = 0.49196 and 
LU2 = 0.09463. The latter have also been confirmed numerically by other codes |5ll . [52l |. We obtain a similar result 
with the code based on Boyer-Lindquist coordinates. 

Quasinormal mode frequencies are the same for all observers. Including null infinity in the computational domain 
does not give us a different picture. The tail decay rates, however, vary from observer to observer as discussed in the 
next subsection. 
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B. Tail decay rates 



111 contrast to the quasinormal mode frequencies, the asymptotic tail decay rates at a finite distance and at null 
infinity are different. Also, before the asymptotic decay rates are obtained, the signals differ substantially from each 
other. The signal that is relevant for an idealized observer is the one measured at null infinity. This is because 
realistic sources of gravitational radiation are at astronomical distances from us. Our detectors are better modeled 
by observers at future null infinity than by observers at finite distances in numerical computations (4ll . [bsI . [66j . 
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Figure 3: Tail decay rates for initial data with a pure / = mode in ingoing Kerr (left) and Boyer-Lindquist (right) coordinates. 
The solid blue curves are measured by observers along J"^ . The solid red curves are measured by observers close to the black 
hole at r = 10 in units where mass m = 1. The dashed black curves correspond to observers at r = {5000, 1000, 500} from top 
to bottom. 



The transition between timelike and null infinity in the decay rates for a monopole perturbation can be seen in 
Fig. [3l In that figure, the local power index has been plotted for various observers ordered by distance ranging from 
close to the black hole up to and including future null infinity. We define the local power index as 



dln|(/)(r,p,y)| 
dlnr 



(22) 



We calculate the time derivatives by using finite differencing from numerical data. The function "Pp.yij) becomes the 
exponent of the polynomial decay of the solution asymptotically in time. We did not observe a significant variation 
of the index with respect to the angular location y of the observer. 
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Figure 4: Tail decay rates for an i = 1 mode. Observers are located as in Fig. |3] 



We see that the monopole perturbation decays with rate —3 near timelike infinity, while the rate at null infinity is 
—2. The decay rate for the I ~ 1 mode is —5 near the black hole and —3 near null infinity, as can be seen in Fig. |4l 
This behavior is the same as in Schwarzschild spacetimc and is in accordance with the prediction of Table II in [55| . 
The decay rate for ^ = 2 is as for I = and the rate for ^ = 3 is as for I ~ 1 due to mode coupling in accordance with 
earlier predictions and calculations [13, [H, [13i UHl ■ For the I = 4 mode the decay rate near timelike infinity is —5 but 
along null infinity it is —4. This is also in accordance with Hod's prediction (see Table I in [ssj). 

We see that the calculation in ingoing Kerr coordinates near the black hole is not as accurate as in Boyer-Lindquist 
coordinates. This is due to the smaller time step in the ingoing Kerr case as mentioned in Sec. IV Al 

We mention that from a physical point of view, in an axisymmetric spacetimc the concept of a pure Z-mode initial 
data does not have a geometrically invariant meaning. In this case, only the monopole perturbations can be defined 
invariantly. This has led to some confusion about decay rates for Kerr spacetimc. In our studies the asymptotic decay 
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rates in both coordinate systems coincide for any /-mode initial data as in ([^D|) due to the fact that we are using the 
same type of initial data in coordinates that are related in a way which leaves the notion of an Z-mode invariant as 
explained in [s^l- Here, we found that the asymptotic rates are the same for such initial data in ingoing Kerr and 
Boycr-Lindquist coordinates also at null infinity. 




500 1000 1500 2000 2500 1000 2000 3000 4000 5000 

Figure 5: Tail decay rates for an I — 4 mode. Observers are located as in Fig. |3] 

Predictions of late-time decay rates are valid asymptotically. In our numerical simulations, we also observe the 
intermediate behavior of the local power index. It is well known that for distant observers the limiting rate is attained 
later in time [HJ. One observes that the decay rate for distant observers is between the rate for close observers and 
the rate along null infinity [4l|, [6^ . This is also seen for the lower Z-modes in Kerr spacetime (Figures |3] and 2]) . For 
the I = 4 mode, however, we observe that the decay for distant observers is faster than for observers very close to the 
black hole (Fig.[5|). We do not have analytical evidence explaining this behavior. It may be su gge sted that it is due to 
the anomalous nature of ? > 4 modes related to mode coupling in Kerr spacetime [10, Issl . Is^. l55l| . We speculate that 
for distant observers the Schwarzschild decay rates may still have some validity. In Schwarzschild spacetime, a decay 
rate of —4 at null infinity would lead to a decay rate of —7 along timelike surfaces. This rate seems to be approached 
by distant observers in Fig. [5l We expect that the rate for close by observers will be attained by distant observers 
much later in time. We wish to stress, however, that this is a speculative attempt at explaining what we observe. 
Analytical calculations need to be performed to understand the intermediate behavior of the local power index. 

VI. DISCUSSION 

We have presented the first numerical evolution scheme for scalar fields on the background of a rotating black 
hole that includes future null infinity. Our method is based on the construction of hyperboloidal surfaces in Kerr 
spacetime presented in [sBl- While previous studies of characteristic foliations do allow, in principle, compactification 
at null infinity on a Kerr spacetime [46l. [47l. [48l. [49| . they did not culminate in numerical applications. In this paper 
we have presented and implemented two methods based on the hyperboloidal compactification technique that include 
null infinity in the computational domain. Our methods are very fiexiblc and allow us to keep any preferred choice of 
coordinates near the black hole. 

We confirmed well-known results concerning quasinormal mode frequencies and tail decay rates measured by ob- 
servers close to the black hole. In addition, we presented decay rates at and near null infinity confirming earlier 
predictions [55^ for the first time. Our method allows us to study the transition between timelike and null infinity 
with respect to the decaying signal. 

We expect that our methods can be applied to discuss physically more interesting evolution equations so that 
predictions concerning the behavior of electromagnetic and gravitational perturbations at future null infinity can be 
studied as well [H, [63|. The coordinate based approach of Sec. Ill may be easier to apply to problems where the 
conformal transformation behavior is not known. In numerical calculations the equations of interest are given anyway 
in a specific coordinate system. Once the asymptotic form of the height function is known, the coordinate based 
approach can be applied easily to transform the equations to include null infinity in the computational domain. 

In numerical applications of the hyperboloidal approach, one usually finds that the hyperboloidal method is very 
efficient in calculating the outgoing radiation accurately [4l|, In our study of spacelike matching, however, we 
find that the presence of a matching domain is a significant source of numerical error. There are many possibilities to 
solve this problem. One option is to devise better matching techniques, such as matching at a multipatch boundary. 
This approach will still have to deal with the change of characteristic speeds across the matching domain. The other 
option would be to construct an analytic, horizon penetrating hyperboloidal foliation which renders the introduction of 
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matching domains unnecessary. It would be interesting to study hyperboloidal foliations of Kerr s pac etime analytically 
or numerically using methods that proved useful in the case of Schwarzschild spacetime [tO, ItH. l72j| . 

This work suggests that once the spherically symmetric case has been mastered, extension to less symmetries is 
more or less straightforward in the hyperboloidal approach. There are two reasons for this. First, hyperboloidal 
foliations can be made to coincide with any known foliation near the black hole [3^. This allows us to use methods 
that have been applied successfully in numerical relativity in an inner domain where the fields arc strong. Second, 
asymptotically flat spacetimes become flat, and therefore spherically symmetric in the asymptotic domain. As a 
consequence the spherically symmetric case already gives important clues on how to deal with the asymptotic region 
in the hyperboloidal approach. 

The observations presented in this paper may be useful in the fully nonlinear case, once the hyperboloidal evolution 
problem with the Einstein equations is solved. Specifically, the numerical discussion of having null infinity at an 
arbitrary but fixed spatial coordinate location and setting the conformal factor to unity in the interior while allowing 
it to vanish in the asymptotic domain in a certain way can be expected to be relevant in the nonlinear case. 
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